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ABSTRACT 

We report results of a Hubble Space Telescope (HST) campaign with the 
Advanced Camera for Surveys to observe Europa at eastern elongation, i.e. Eu- 
ropa's leading side, on 2008 June 29. With five consecutive HST orbits, we 
constrain Europa's atmospheric O I 1304 A and O I 1356 A emissions using the 
prism PR130L. The total emissions of both oxygen multiplets range between 
132 ± 14 and 226 ± 14 Rayleigh. An additional systematic error with values 
on the same order as the statistical errors may be due to uncertainties in mod- 
elling the reflected light from Europa's surface. The total emission also shows a 
clear dependence of Europa's position with respect to Jupiter's magnetospheric 
plasma sheet. We derive a lower limit for the O2 column density of 6 x 10 18 
m~ 2 . Previous observations of Europa's atmosphere with STIS in 1999 of Eu- 
ropa's trailing side show an enigmatic surplus of radiation on the anti- Jovian side 
within the disk of Europa. With emission from a radially symmetric atmosphere 
as a reference, we searched for an anti- Jovian vs sub- Jovian asymmetry with 
respect to the central meridian on the leading side and found none. Likewise, 
we searched for departures from a radially symmetric atmospheric emission and 
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found an emission surplus centered around 90° west longitude, for which plausi- 
ble mechanisms exist. Previous work about the possibility of plumes on Europa 
due to tidally-driven shear heating found longitudes with strongest local strain 
rates which might be consistent with the longitudes of maximum UV emissions. 
Alternatively, asymmetries in Europa's UV emission can also be caused by in- 
homogeneous surface properties, inhomogeneous solar illuminations, and/or by 
Europa's complex plasma interaction with Jupiter's magnetosphere. 

Subject headings: Europa, Jupiter, atmospheric emission 
Version: January 27, 2013 



1. Introduction 



Jupiter's satellite Europa is an outstandingly interesting body in our solar system as 
several lines of evidence point to a subsurface water ocean under its icy cr ust. Europa is a 
differentiated body with an outer layer of H2O with a thickness of ~100 km ([Anderson et al. 



20071 ; iHussmann et al.ll2002l ). Its surface is relatively young with an average age of ~60 Myr 
and it is much less heavily cratered than its two neighboring satellites Ganymede and Callisto 
(e.g., iGreeley et al.l 120041 ) . Evidence for Europa's subsurface ocean com es from observations 



of geological structures (e.g. , iPappalardo et al.lll999t IGreeley et al.ll2004[ ) . Europa's surface is 



rich in various features, such as ridg es, chaotic regions, and cycloids possibly due to fl exing of 
its icy crust under tidal forces (e.g., IPappalardo et al.lll999l ; iNimmo fc Gaidosl 120021 ) . These 
features are naturally explained by a subsurface ocean of liquid water, but it cannot be ruled 
out that they are alternatively du e to processes in warm, soft ice with only localized or partial 
melting ( IPappalardo et al.lll999[ ). A strong and fully indepe ndent argument for a subsurface 



ocean comes from observations of induced magnetic fields (IKhurana et al. 



1998; 


Kivelson et al. 


2000; 


Zimmer et al. 


2000; 



1998 



201 



Neubauer 
oj) which 



require a subsurface layer with a sufficiently large electrical conductivity Within Europa's 
geological and geophysical context, only saline liquid water is considered a viable candidate 
to achieve the required conductivity. 

Europa also possesses a thin oxygen atmosphere discovered w ith the Goddard High- 
Resolution Spectrograph (GHRS) of the Hubble Space Telescope by lHall et al.l ( 119951 ). Ob- 



served fluxes of 11.1 ± 1.7 x 10 5 photons cm" - 2 s- 1 and 7.1 ± 1.3 x 10" 5 photons cm 
s _1 (at the Telescope), corresponding to a brightness of 69 ± 13 and 37 ± 15 Rayleigh, for 
the two oxygen multiplets O I 1356 A and O I 1304 A , respectively, im plied a very thin 
molecular oxygen atmosphere with a surface pressure of about 10~ n bar (jHall et al.l 11995 
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19981 ). The radiation is excited predominately by electron impact dissociation of molecular 
oxygen. The previous observations were taken on 1994 June 02 at western elongation, i.e. on 
Europa's trailing hemisphere. Subsequent observations on 1996 July 21 with GHRS on the 
leading hemisphere rendered fluxes of 8.6 ± 1.2 x 10~ 5 photons cm -2 s' 1 (54.7 ± 7.6 R) for 
O I 1304 and 12.9 ± 1.3 x 10' 5 photons cm' 2 s' 1 (82.1 ± 8.3 R) for O I 1356. Observations 
on 1996 July 30 confirmed the oxygen atmosphere on the trailing hemisphere with 12.6 ± 
13 x 10~ 5 (80.2 ± 8.3 R) for O I 1304 and 14.2 ± 1.5 x 10~ 5 (90.4 ±9.5 R) for O I 1356 
( lHall et al.lll9981 ). The joint fluxes of both oxygen multiplets were thus 18.2 x 10~ 5 photons 
cm -2 s' 1 (106 R) on the trailing side in 1994, 26.8 x 10' 5 photons cm' 2 s' 1 (170 R) also 
on the trailing side in 1996, and 21.5 x 10~ 5 photons cm ' 2 s' 1 (137 R) on the leading side 



in 1996. Based on these observations, lHall et al.l ( 119981 ) estimated O2 column densities in 



the rang e of 2.4 to 14 x 10 m . C omprehensive reviews on Europa's atmosphere can be 



found in iMcGrath etaP tool 12009^ 



Europa's atmosphere is considered to be primarily generated through surface sputter 



ing by magnetospheric ions ( 



ohnson et al.l Il982l ; iPospieszalska fc Johnson! Il989l ; IShi et al. 



1995c IShematovich et al.l 120051 ) . Dominant loss processes are atmospheric sputtering and 
ionization. The balance of production and loss rates, which both also depend on the elec- 
tromagnetic field environment establishes the average atmospheric content. Estimations of 



this balance leads to a O2 c olumn density of ~5 x 10 18 m 2 and a neutral O2 loss rate of 



50 kg s 1 (jSaur et al.lll998l ). The neutral loss, includin g H? from thermal e scap e, generates 



a neut ral torus along the orbit of Europa discovered by iMauk et al.l (120031 ) and lLagg et al. 
(|2003h . 



Europa is embedded in the Jovian magnetosphere whose plasma constantly streams 
past Europa and its atmosphere. The interaction of the magnetospheric plasma with Eu- 
ropa's atmosphere is primarily responsible for the generation of an ionosphere on Europa 



dKliore et al. 



19971 ) . This interaction i n addition to induction in a possible sub-surface ocean 



creat es magnetic field perturbatio ns dKivelson et a 



20081 ) and generates Alfven wings (INeubauer 



199C 



1997 



Saur 



Saur et al.lll998l ; ISchilling et al. 



2004]) , which lead to auroral foot - 



prints of Europa in Jupiter's upper atmosphere (IClarke et al.l 120021 ; iGrodent et al.l 120061 ). 
Europa's plasma interaction modifies not only the electron density but also the electron 
temperature in Europa's atmosphere. Both the electron density and temperature control 
the O 11304 A and O I 1356 A e mission rates. Europa's plasma inter act ion has been studied 
by ISaur et all f)l998h : iLiu et al.l (boOlh : ISchilling et all koQ'A l2008h . where the latter two 
models self-consistently include induction in Europa's possible subsurface ocean. 



After the initial spectral observation of Europa's atmosphere by lHall et al.l ( 119951 . 119981 ). 
the first spatially resolved observations of Europa's atmospheric emission were performed 
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with the Space Teles cope Imaging Spectrograph (STIS) in 1999 October 5 during five con- 
secutive HST orbits (IMcGrath et al.ll2004 I2009T ) . A superposition of all exposures is shown 
in Figured! The emission displays a clear asymmetry with respect to the sub- Jovian/anti- 
Jovian side. On the anti- Jovian side within the disk of Europa occurs a significant enhance- 
ment of the emission. Even though Figure [TJ suggests time-stationarity of the emission, 
the location of the maximum of the emission moves from Europa's norther n hemisphere to 



its so uthern hemisphere, but stays predominantly on the anti- Jovian side (IMcGrath et al. 



20091 ). This emission is puzzling as it occurs within the disk of Europa, while for a radially 
symmetric atmosphere combined with a radially symmetric plasma interaction, t he emission 
would be largest near the limb just outside of Europa's disk. ICassidy et al.l (120071 ) investigate 
such a stationary surplus of the oxygen emission and show that a non-uniform distribution 
of reactive species in Europa's porous regolith can result in a non-uniform O2 atmosphere 
consistent with Figure [TJ 

In this work we investigate HST/ ACS observations of Europa's UV emission taken on 
2008 June 6 when Europa was at eastern elongation to further constrain Europa's atmo- 
spheric content and to investigate if there are also enigmatic inhomogeneities in Europa's 
UV emission on the leading hemisphere. In section [2] we present our observations at eastern 
elongation and describe how they were analyzed. In section [3] we present the results of our 
analysis and discuss their implications. Finally, in section HJ we summarize our results and 
point to unresolved issues. 



2. Observations and Analysis 



During HST program 11186, five consecutive orbits of observations with ACS/SBC 
prism PR130L on 2008 June 29 were taken (see Table [TJ). Europa was at eastern elongation, 
i.e. the leading side of Europa was visible from Earth. Each orbit was subdivided into four 
distinct exposures of ~10 min. During the time of the observations Europa also crossed the 
current sheet of Jupiter's magnetosphere. Europa's positions with respect to the current 
sheet are indicated in Table [TJ 

Due to th e unavailability of STIS during the time of our observations, we used the 
prism PR130L ( iBohlin et al.ll2000l ; lLarsonll2006l ) to separate reflected solar photons from the 
surface of Europa and the oxygen multiplets O I 1356 A and O I 1304 A from Europa's 
atmosphere. The prism covers nominally the wavelength range larger than 1216 A to ~2000 
A and thus suppresses Ly a. Unfortunately, the wavelength range larger tha n 2000 A is less 
stron gly suppressed than anticipated leading to a significant 'red leak' (e.g. iFeldman et al. 



201 ll ). The sensitivity and dispersion properties have been recalibrated on-orbit by 



Larson 
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( 120061 ). The dispersion, i.e. A A per pixel, is a highly non-linear function of wavelength. The 
wavelength A as a function of trace distance Ax on the detector can be written in the form 



A(Ax) = a x + a 2 (Ax — a ) 1 + a 3 (Ax — a ) 2 + a 4 (Ax — a ) 3 + a 5 (Ax 



with the coefficients a t o a?, depending on t he object po sitions x re f,y re f in the direct im- 
age, i.e.. without prism (jBohlin et al.ll2000t lLarsonll2006l ). The spectral trace is addition- 



Bohlin et al. 


(2000 


) and 


Larson 


(2006) 



In Figure [2] we display the sum of the four exposures in each orbit separately. A su- 
perposition of all exposures is shown in the top panel of Figure SJ In each image, celestial 
north is toward the right and defines the x-direction considered in the present study. The 
top of the images is towards the anti- Jovian side of Europa and is labeled y-direction. The 
dispersion of the prism is in the x-direction, i.e. longer wavelengths lie further to the right 
in each image. In each image in Figure [2] and HI the two ellipses show the expected locations 
of Europa's disk for the O I 1304 A and O I 1356 A emissions observed with the PR130L 
prism. Unfortunately, both images overlap. The overlap is even somewhat larger than in- 
dicated in the overlapping ellipses as the radiation from Europa's atmosphere off the disk 
additionally contributes to the radiation. Therefore we analyze the sum of the O I 1304 A 
and O I 1356 A emission in this paper and investigate possible asymmetries of the joint 
emission in the y-direction only, i.e. the sub- Jovian/anti- Jovian direction, as this direction 
is not affected by the overlap. 

The sunlight reflected by Europa's surface significantly contributes to the measured 
fluxes. In particular, emission longward of ~2000 A is not fully filtered out and leads to 
the triangle shaped emission on the right hand side of the image displayed in Figure [2J This 
long wavelength 'red leak' exceeds the fluxes of the O I images by roughly two orders of 
magnitude much further to the right of the O I images as displayed in Figure [5j 

In order to account for the contribution of the reflected sunlight to the O I images, we 
model its photon flux F s as a function of position (x, y) on the detector. The reflected and 
dispersed solar light can be calculated as 



F s (x, y)= dX p(X)t(X)a(X)D[x - (Ax(A) +x E ),y- (y d (Ax(X)) + y E )\ (2) 

./solar spect 

with the dispersion dX/d(Ax) ( Larson! bood ) and Ax(A) being the trace distance of the 
dispersed photon flux in x-direction as a function of wavelength A compared to the location of 
the undispersed flux without prism . The location of the spectral trace is also shifted in the y- 
direction by yd{Ax) (jLarsonll2006l ). The coordinates x E , and y E refer to the center of Europa's 
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disk without dispersion. D[x(X), y(X)] represents the disk of Europa with D[x, y] — 1 within 
the disk, and D[x,y] = outside of the disk. Since we used a pr i sm, D is a function of 
wavelength A. The solar photon flux p(A) is given by I Woods et al.l ( 119961 ) and displayed in 
Figure El t(A) is the throughpu t, and a(A) Europa's albedo. The throughput for A < 2000 
A is provided by |L arson! ( 2006 ) and constrained for A > 2000 A based on observations of 16 
CygB as detailed in Appendix IA.1I The flux F s (x,y) is additionally convolved with a point 
spread function as described in Appendix IA.2I 

Since no direct image of Europa was taken during the HST campaign 11186, the location 
of the direct image has to be reconstructed from the dispersed images. For this purpose the 
strong red leak is very useful. The long wavelengths reflected solar light generates in the 
obtained observations an easily discernible region of maximum flux corresponding to the size 
of Europa when plotted with a linear color scale, owing to the very low dispersion of the 
prism at red leak wavelengths. For this long wavelength image of Europa, we determine 
its coordinates and then use expression fl2]) to subsequently determine xe and yE- The 
latter step is possible despite the uncertainties in the throughput for long wavelengths since 
the dispersion is very large at long wavelengths, e.g. ~ 250 A pxl" 1 at 4000 A. Therefore 
uncertainties in the slope of the throughput at long wavelengths have very small effects on 
the position of the dispersed long wavelength image of Europa. 

Europa's albedo as a func tion of wavelength is only partially constrained by previous 
observations. lHall et al.l ( 119981 ) derive a disk averaged albedo of 0.015 ± 0.004 on the leading 
side and similar albedo values on the trailing side at A = 1335 A . Based on an analysis of 
STIS observations of Europa's trailing side in 1999 (ID 8224, PI M. McGrath), the albedo 
stays remarkably flat in the FUV between 1400 and 1700 A at values near 0.016. The 
constancy of the albedo values in this wavelength range (not shown here) is d erived similarly 



to the an alysis of STIS obser vations of Ganymede (ID 7939, PI W. Moos, iFeldman et al. 



( 120001 )). iHendrix et al.l ( 119981 ) find that the albedo on the leading side at central meridian of 
12 0° W increases ap prox imately linearly from ~0.2 to ~0.65 within A = 2200 A to A = 3200 
A. INoII et al.l ( 1995 ) and IHendrix et al.l ( 1998 ) find a lower albedo on the trailing side which 
changes from ~0.08 to ~0.17 within A = 2200 A to A = 3200 A. The general t rend that the 



albedo on the leading side is larger compared to the trailing side is also seen by iNelson et al. 
( 119871 ). The authors observe albedo variations by roughly a factor of two from the leading 
to the trailing hemisphere within A = 2400 A to 3200 A. The maximum values of the albedo 
on the leadin g side are 0.18 (2400-2700 A), 0.24 (2800-2900 A) and 0.32 (3000-3200 A). 
Hansen et al.l ( 120051 ) derive an albedo of ~1% (without error bars) near 1335 A and at 94° 
phase angle. 



Our analysis pertains to the leading side, i.e. for the central meridian 90°, where we 
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assume the model albedo shown in Figure El In the absence of albedo measurements within 
1700 A and 2200 A, we choose to construct an albedo model which is as simple as possible 
and is still consistent with the constraints described in the previous paragraph. For A < 1700 
A our model albedo is constant at a value of 0.015. For 1700 A < A <3000 A the albedo 
rises linearly from 0.15 to 0.41 and stays flat for A >3000 A (see Figure [6]). Further details 
on how the quantitative values of the linear rise of our model albedo are chosen, will be 
discussed in subsequent paragraphs where we compare the observed fluxes with the modeled 
solar reflected light. 

With the procedure and model inputs described in the previous paragraphs, we calculate 
how the solar light reflected from the surface is dispersed with the prism PR130L and appears 
on the detector. The resultant flux is shown in Figure IU lower panel. For comparison we 
show the observed fluxes of all exposures superposed in Figure HI top panel. Several features 
of the modeled reflected light can be seen in the observations, such as the bright emission 
'triangle' to the right of the O I emissions due to reflected long wavelengths light. The 
extension of the emission in y-direction within the triangle is due to the scattering within 
the prism and described by the point spread function constrained in Appendix IA.2I It is 
also clearly visible that there is a contribution from the reflected light also within the O I 
images. 

In Figure [5] we quantitatively show the modeled reflected surface flux and the observed 
flux of all exposures superposed along the direction of the dispersion. To reduce statistical 
fluctuations, we integrated the fluxes in the y-direction. To minimize the contribution from 
the off-disk emission, the integration was performed over 25 pixels from rows 361 to 385 (i.e., 
somewhat less than the diameter of Europa on the detector, which is approx. 30 pixels in 
the y-direction). The integrated modeled (black) and observed fluxes (red) are plotted as a 
function of trace distance (i.e., in the x-direction) in Figure [5J The two vertical dashed lines 
show the locations of the O I 1356 A image and the two vertical dotted lines the locations 
of the O I 1304 A image corresponding to the size of Europa's disk. The solid vertical line 
in Figure |5] shows the 2000 A 'influence line'. Reflected light from the disk of Europa with 
wavelengths larger than 2000 A, i.e., from the red leak, would lie strictly to the right of this 
line if scattering were not to take place. 

We use the observed fluxes displayed in Figure [5] to constrain our model albedo (see 
Figure [6]). We assume a simple linear increase of the model albedo for wavelengths larger 
than 1700 A. The quantitative values of this increase has been chosen to match the observed 
fluxes sufficiently far away from the wavelength ranges which correspond to the O I emission 
from Europa's atmosphere (i.e. outside of the vertical dotted and dashed lines). The modeled 
reflected light is shown as a solid black line in Figure which fits the observations outside of 
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the atmospheric O I wavelengths ranges in general well. However, near the solid vertical line 
the surface reflected light overestimates the observed fluxes and does not reproduce the 'knee' 
in the observations. A principal reason could be a dip in the albedo around 2000 A. But 
any uncertainty in the model albedo a can equally well be reflected in an equal uncertainty 
in the throughput t when working with the prism observations. Only the product of a and 
t constrains the observed fluxes (see equation (J2])). The throughput at wavelengths around 
and larger than 2000 A is also not well constrained (see discussion in Appendix I A . 1 [) . To 
improve the fit of the modeled reflected light we introduce a throughput dip of a factor 
1/3 within 1800 A and 2300 A (shown as dashed blue line in Figure [T5|) . This modified 
throughput leads to a modified spectrum of reflected and dispersed light from the surface 
of Europa and is displayed as the dashed black line in Figure [5j The spectrum with the 
throughput dip reproduces the knee in the observed spectrum much better. We will call it 
in the remainder of this paper 'alternative' model of reflected surface light. It will provide 
in subsequent parts of the paper a rough estimate for the consequences of the uncertainties 
in the model albedo and/or throughput, i.e., it provides an estimate for systematic errors 
introduced in the subtraction of the reflected light. 

At some distance away from the locations of the oxygen images (constrained by the 
vertical dashed and dotted lines) the observation are close to our modeled reflected light. 
The surplus of the observed emission compared to the model solar reflected light within and 
near the dotted and dashed lines is the O I 1356 A and O I 1304 A radiation emitted from 
Europa's atmosphere. It should be noted that in Figure the flux scale is logarithmic, thus 
a surplus at x=280 appears roughly ten times larger than the same surplus at x=330. 

We note that any possible contribution by the geocorona plays a negligible role in our 
analysis. The O I 1356 A of Earth's atmosphere is below the HST orbit. As we observe 
Europa near Jupiter opposition, i.e., during HST orbit night, the O I 1304 A emission only 
contributes 1 or 2 Rayleighs uniformly on the detector. 

In the next step we display in Figure [7J the observed fluxes for each HST orbit separately. 
For comparison we also show the modeled solar reflected fluxes. The fluxes are integrated 
between rows 361 to 385 similar to Figure EJ All error bars include statistical errors only. 
The signal to noise ratio is calculated by S/N = (S — B — D)/yS + B + D with S and 
B corresponding to the total counts of the signal and the background, respectively, and D 
to the total counts of the ACS/SBC dark noise (S, B and D are summed over exposure 
time and pixel u nder consid eration, respectively). The dark noise count rate is 1.0 xl0~ 5 



cts -1 s _1 pixel -1 ( ICoxl 120041 ) . A more detailed representation of the emission from Europa's 
atmosphere is given by the integrated fluxes on a linear scale spanning wavelengths below 
2000 A. The solid vertical line in Figure [7] shows the 2000 A 'influence line'. The two vertical 
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dashed lines show the locations of the O I 1356 A image and the two vertical dotted lines 
the locations of the O I 1304 A image corresponding to the size of Europa's disk. 

The differences between the modeled surface reflected light (black) and the observed 
fluxes (red) in Figure [7| are due to the emission from Europa's oxygen atmosphere (blue). 
We multiplied the values of emission from Europa's atmosphere by a factor of five to more 
clearly display it within the individual plots. The corresponding values for the atmospheric 
emission are shown on the right hand side of each panel in blue. The model fluxes (black) 
to the right of the O I emission region and to the left of the 2000 A influence line match all 
individual exposures and the superposed exposure (red) fairly well and thus demonstrate that 
the solar reflected light contributions does not change from orbit to orbit. In contrast, the 
observed emissions in blue display apparent time-dependence. The observed fluxes change 
as a function of orbit number as will be analysed in more detail in the following section. The 
observations also show that there is excess emission beyond the disk of Europa due to the 
spatial extension of Europa's atmosphere as visible to the left of the O I 1304 A disk, i.e. to 
the left of the left dotted vertical line. Any excess emission to the right of the O I 1356 A 
disk, i.e. to right of the right vertical dashed line, might be buried within the large error bars 
in this region. The signal to noise ratio of the radiation from Europa's atmosphere displayed 
in Figure [7| is however partially as low as two as apparent in the large error bars. Therefore 
we will investigate in the remainder of this work areas corresponding to a larger number of 
pixels for a better signal to noise ratio. 



3. Results 

3.1. Radial/equatorial structure of atmosphere 

In Figure [S] we show the average emission as a function of radial distance from the center 
of the disk of Europa, which we calculate as average emission within thin radial shells of 
both disks (O I 1304 A and 1356 A). The emission is normalized to the number of pixels in 
each shell. The radiation within the shells is calculated as the difference in radiation within 
radius r + dr and of radius r. For radial values where the two O I 1304 A and 1356 A areas 
overlap, we do not calculate the overlapping part twice. No attempt is made to separate the 
O 1 1304 emissions from the O 1 1356 emissions in the overlapping regions. Instead our model 
instrument response replicates the overlap as part of the fitting process. We restrict ourselves 
to pixels within columns 274 to 323 in case the shells extend beyond these columns. These 
columns just embrace the O I 1304 A disk to the left and the O I 1356 A disk to the right 
in x-direction. If we were to integrate further to the right in the x-direction (direction of the 
dispersion) we would also integrate pixels whose emission is significantly contaminated by 
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the red leak. To avoid over-counting the O I 1304 A compared to the O I 1356 A emission, 
we do not integrate further to the left of the O I 1304 A image. Consequently, the profiles 
in Figure [H] is more constrained by equatorial than by polar off-disk emission. 

The red line in Figure [8] shows the original observations. The dotted line shows the 
reflected surface light without inclusion of the point spread function. The apparent model 
emission extends to radial distances larger than 1 R E even though the surface reflected light 
is generated within 1 R E . The reason is that the dispersion of the continuous solar spectrum 
also generates emission out side of the location of the O I images. The solid black line 
shows the model surface reflected light including the point spread function. It shows that 
the scattering of photons within the prism is a significant effect as can be seen in Figure H] as 
well. The model reflected light fits the observations at large radial distances remarkably well 
and thus independently confirms the quality of the constructed point spread function (see 
Appendix). The difference between the solid red and the black line is the excess emission from 
Europa's atmosphere (shown as dashed solid line). It demonstrates non- negligible emission 
above the disk of Europa which tapers off towards ~1.5 - 1.7 Rj. In Figure [8] we also show the 
resultant flux from Europa's atmosphere calculated with the alternative throughput model 
as red dashed lines. These fluxes differ by only ~10% or less from the fluxes derived with 
our standard model. 



3.2. Atmospheric emission, individual orbits 



The O I 1304 A and the O I 1356 A emission from Europa's atmosphere can be re- 
constructed by subtraction of the reflected and dispersed solar light of the surface from the 
observations (see Section [2]). The emissions from Europa's atmosphere for all five orbits and 
the sum of all exposures are shown in Figure [9j The color scale is similar to Figure [2], but 
enhanced by a factor of 1.5 to more clearly display the flux variations. We note that individ- 
ual pixels or small cluster of pixels should be interpreted with much caution due to their low 
signal to noise ratio. The general impression from the images in Figure M is that most of the 
flux occurs within or close to both O I images. The images also show time dependence with 
orbit 2 and 3 having maximum fluxes. The emission emerges mostly within the disk of Eu- 
ropa, but it also extends b eyond the disk as due to a finite atmos pheric scale heig ht reported 



in pr evious observations (jHall et al.l Il995l ; iMcGrath et al.l 120041 ) and modeling (ISaur et al. 
1998h . 



- 11 - 



3.3. Total fluxes 

Due to the low signal to noise ratio for individual pixels or small groups of pixels, it is 
necessary to integrate the emissions over larger areas. Therefore we calculate the total flux 
from Europa's atmosphere for each orbit within an extended disk including a radial distance 
of 1.3 Re from Europa's center with Europa's radius Re = 1569 km. In Figure [TTJ1 the total 
fluxes are shown in red as a function of Europa's position with respect to the current sheet 
of Jupiter's magnetosphere. Jupiter's magnetic moment is tilted with respect to Jupiter's 
spin axis by ~9° . Jupiter's fast rotation period of 9 h and 55 min is responsible for large 
centrifugal forces, which confine the magnetospheric plasma into a disk around Jupiter's 
centrifugal equator, called 'current sheet' or 'plasma sheet'. The current sheet is close to the 
magnetic equator while Europa's orbital plane is within Jupiter's spin equator. Thus the 
current sheet rocks above and below Europa within 11.23 h (the synodic rotation period of 
Jupiter as seen from Europa). Therefore Europa is exposed to different plasma conditions, 
in particular varying electron densities, within the 11.23 h. This effect is clearly visible in 
Figure [TU1 The maximum emission of ~36 ± 2.2 x 10 -5 photons cm' 2 s _1 is reached when 
Europa is in the center of the current sheet and minimum values of ~21 ±2.1 x 10~ 5 photons 
cm -2 s _1 occur when Europa is at the edge of the current sheet. These fluxes correspond 
to ~230 ± 14 and ~130 ± 14 Rayleigh, respectively (see also Table 1). The quoted error 
bars contain statistical errors only. In Figure dUJ we also show the fluxes calculated with 
the hypothetical dip in the throughput curve as alternative model with dashed lines. The 
alternative fluxes differ from the standard model by as much as the width of the statistical 
errors. The difference between the standard and the alternative model can be considered 
as a rough measure for the systematic error introduced by subtracting the surface reflected 
light, which contains uncertainties due to the partially unconstrained albedo and red leak. 

In black we show the fluxes within the disk of Europa, i.e., within 1 Re- The associated 
values are roughly 3/4 of the fluxes of the extended disk. We note, as both O I images 
overlap, some of the off disk emission from the real atmosphere falls within the O I disks 
of the two dispersed O I images. Thus the quoted ratio is an upper limit for the real disk 
emission. 

We can compare the current sheet dependence of the tot al fluxes to the observation s 



from 1999 October 5 taken with STIS on the trailing side flMcGrath et all I2004L 120091 ) . 
Similarly to the ACS data, we also integrate the STIS data within 1.3 Re- The total fluxes 
from the 1999 October 5 campaign are shown as brown diamonds in Figure [10J The green 
squares in Figure [10] show the current sheet dependence of the O I 1356 A emission and 
the violet triangles those of the O I 1304 A emission. The ratio of the latter two is roughly 
two to one. The 1999 observations also show a clear dependence on Europa's position with 
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respect to the current sheet for the O I 1356 A emission, while the dependence of the O I 
1304 A emission is ambiguous. Their total fluxes are roughly 30% to 40% larger than the 
values derived from the 2008 ACS observations of this work. The combined O I 1356 and 



O I 1304 fluxes of 137 Rayleigh observed by lHall et all fll998f ) with GHRS on the leading 
side are consistent with the 2008 ACS observations outside of the current sheet and ~50% 
smaller than the 2008 ACS values when Europa is in the center of the current sheet. 

We can use the observed fluxes within the range of 170 to 270 Rayleigh and estimate 
average column densities. Assuming for simplic it y a c onstant electron density of 40 cm -3 
and temperature of 20 eV similar to lHall et al.l ( 119951 ) and neglecting the electrodynamic 
interaction, we can derive a lower limit for the O2 column density of 6 to 10 x 10 18 m -2 . 
Our derived column densities show somewhat less variations, but are still roughly co n sisten t 



compared to the column densities of 3.5 to 11 x 10 m derived by lHall et al.l ( 119981 ). 



Note that the observ ed electron density upstream of Europa varies between 18 to 250 cm 3 
kivelson et~aDl2004h . 



3.4. Search for asymmetries on leading hemisphere 



The spatially resolved observations of Europa's atmosphere on the orbital trailing hemi- 
sphere in 1999 show a pronounced surplus of emission on the anti- Jupiter side (see Figured]). 
We quantitatively searched for similar asymmetries on Europa's leading hemisphere. Fortu- 
nately, the dispersion is in the north-south direction and we can thus well analyze possible 
asymmetries with respect to the sub- Jovian and anti- Jovian side. In Figure [TTJ we show 
the ratio of the total flux within the sub- Jovian side compared to the anti- Jovian side of 
the emission for each orbit. The black triangles are calculated from the emission within the 
two O I disks only, while the red stars show the ratio calculated from an extended disk, i.e. 
within a radius of 1.3 Re- The sub- Jovian/anti- Jovian emissions within the disk are within 
the error bars consistent with a ratio of one, i.e. consistent with a symmetric emission. The 
extended disk emission show s an average ratio of ~0.9 ±0.2. Similar to UV observations of 
Io's atmospheric emission by iRoesler et al.l (Il999l). one might expect the anti- Jovian side of 
Europa to be brighter due to the Hall effect ( iSaur et al.lll999l . 120001 ) . However, this effect 
is not resolvable within the current error bars. In summary, no sub- Jovian and anti- Jovian 
asymmetry on the leading side of Europa is visible in contrast to its trailing side. 

In a next step we quantitatively search if there is substructure in the emission on 
the disk of the leading hemisphere. Such a substructure might be generated by in ho- 
mogeneities of Europ a's surface properties such as discussed by ICassidy et al.l (120071 ) or 
plumes speculated by iNimmo et al.l (j2007al ) similar to those observed on Enceladus (e.g., 
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2006 



Burger et al. 



2007 



Saur et al. 



Waite et aU200a ISpencer et aU2006l:|Porco et al.ll2006t iDougherty et al.ll2006l ; lHansen et al. 



20081 ) . Therefore we integrate the fluxes within the disk 



parallel to the dispersion, i.e. along the x-axis. The resultant integrated fluxes are shown 
as a function of the sub- Jovian/anti- Jovian direction as red curves plus error bars in thin 
black for all individual orbits and the sum of all exposures in Figure [12] For comparison we 
show as black curves the expected emission s from within the disk for a radia lly symmetric 
atmosphere with a scale height of 145 km (IHall et al.l Il995l ; ISaur et al.l Il998l ) . The effects 
of the point spread function on the radially symmetric model atmosphere is included. The 
model curve for a radially symmetric atmosphere is arbitrarily normalized and displayed 
to demonstrate possible deviations from a radially symmetric atmosphere in the observed 
emission in the sub- Jovian/anti- Jovian direction. The observations for the individual orbits 
show temporal variations in the overall flux. As described in Section 13.31 overall temporal 
variations of the total intensity are consistent with Europa's varying position with respect 
to the current sheet. 

We use the sum of all exposures, which combines the O I 1356 A and the O I 1304 
A emission, to investigate possible spatial substructure within the emission. There is more 
emission near the central meridian (90 degree at eastern elongation) compared to the emis- 
sion from near the limb as the emission near the limb (near pixel 360 and 390) is smaller 
compared to the emission profile expected from a radially symmetric atmospheric emission. 
The observed emission could be consistent with a density enhancement around the 90 degree 
meridian. Only as an example to characterize the amount of gas, an atmospheric substructure 
with a radius of ~250 km and a density enhancement by a factor of two to three compared 
to the background atmosphere would create the observed deviations in Figure [12] compared 
to a radially symmetric atmosphere. Such den sity variations could be consistent with model 
calculations reported by ICassidy et al.l (120071 ) . In the right lower panel in Figure [121 which 
covers the sum of all exposures, we also show the alternative model derived from an alter- 
native throughput curve as a jagged black line close to the jagged red line. It demonstrates 
that the applied modifications in the throughput does not have qualitative implications on 
our derived conclusions. We also note that in Figure [12] lower right panel at least four data 
points close to the limb are within the error bars below the radially symmetric model. We 
could thus also combine two or all four to a joint point with even smaller error bars. 



Nimmo et al.l (l2007bl ) suggested that the vapor plumes discovered at Enceladus were 



caused by tidally- driven shear heati ng, with the vapor produ ction rate dependent on the local 
strain rate. If similar shear heating ( jNimmo fc Gaidosll2002l ) and vapor production occurs at 
Europa, then the local rate of vapor production should depend on the local strain rate. For a 
thin ice shell, the time average of the squared tidal strain rate varie s spatially, with maxima 
at the poles and equatorial minima at 0° and 180° west longitude ( lOjakangas k, Stevenson 
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19891 ). Figure [131 plots the latitudinal average of this quantity as a function of longitude, 
demonstrating that strain rates (and hence predicted vapor production) are largest around 
90° and 270° west longitude. This analysis is probably over-simplified, in that it assumes 
the orientation and areal density of shear-heating features is uniform. To take the spatial 
variation of shear- heating features into account, we also calculated the mean resolved shear 
stress (a proxy for shear heating rate) on ma pped lineaments. Th e stress calculations are 
described in the Supplement ary Information oflNimmo et all ( 2007bl ) and the lineament map 
was based on Figure 15b of iKattenhorn fc Hurfordl (120091 ) . Using this approach, we predict 
a peak in vapor produc tion at 230° west, in st riking agreement with the observations on the 
trailing side reported by lMcGrath et all ( 120041 ) and shown in Figure [T]of this paper. However, 
the loca tion of this peak may be an artifact due to the incomplete imaging c overage of polar 



regio ns (Doggct t et alll2009l . Figure 1.). If the observed peak at 230° west (IMcGrath et al 



20041 . 120091 ) is not a transient feature, and arises because of shear heating, then we predict 
that future imaging will reveal a paucity of polar ridges around 270° west. 

The observed surpluses of O I emission near 90° and in between 200° and 250° west 
longitude need further investigations. For example, if the surplus O I emission were confined 
to 1304 A and not 1356 A, it could indicate an optically thick atmosphere in O atoms, 
which would resonantly scatte r solar 1304 A photons with enhanced emission at the center 
of the disk. lHall et all dl995h used the '2:1 ratio' on O I 1356/0 I 1304 intensity to put 
an upper limit on the O at om column densi ty of < 2 x 10 18 m -2 . But this ratio has been 
observed to be < 2:1 (cf. lHall et all ( 119981 )) and thus solar 1304 A resonance scattering 
by O atoms may be non-negligible. The observed surpluses also might be due to varying 
surface properties, such as porosity or albedo. The electrodynamic interaction of Europa's 
atmosphere with Jupiter's magnetosphere will also generate asymmetries in the interaction, 
however it generally acts to enhance emission on the flanks in contrast to the enhanced 
emission in the wake. These effects warrant further investigations which are however beyond 
the scope of this observational paper. 



4. Summary and Discussion 



Europa's O I 1304 A and 1356 A emissions on the leading side observed on 2008 June 
29 with ACS/SBC prism PR130L show total fluxes in the range of 130 to 230 Rayleigh. 
These values are generally comparable with previous observations. They are slightly smaller 
than fluxes obtained with H ST/STIS on the trailing side in the range of 210 to 370 Rayleigh 
(IMcGrath et all I2004J . 120091 ) , but they are in the range or somewhat lar ger than previous 
GHRS observation of 137 Rayleigh on the leading side (lHall et al.l[l998h . The ACS/SBC 
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and the HST/STIS observations demonstrate a clear dependence of the emission depending 
of Europa's position relative to Jupiter's current sheet. We do not find an asymmetry of 
the atmospheric emission with respect to the sub- Jovian/anti- Jovian side, but find a surplus 
of emission near 90° west longitude compared to a radially symmetric atmosphere. The 
previous STIS observations are consisten t with a surplus of emission on the trailing side in 
the range of 200° to 250° west longitude ( IMcGrath et all 12004 120091 ). A surplus of emission 
could be due to a surplus of atmosph eric gas, which could be consistent with predictions of 
possible plume locations on Europa by lNimmo et al.l (j2007al ) near the same longitudes. Other 
reasons which could contribute to the observed longitudinal enhancement in the UV emission 
might be intrinsic or solar illumination driven spatial variations of Europa's ice prop erties 



controlling its albedo, porosity and sputtering properties (e.g., iCassidy et al.l 120071 ). For 
example, the albedo is known to vary with longitude on Ganymede fo r wavelengths larger 
than 2200 A jNelson et al.lll987l : INoII et al.l[l995[ lliendrix et aDll998h however it appears 



to be constant within the error bars in the FUV at 1335 A (IHall et al.lll9981 ). The emission 



patterns of Europa's atmosphere are also controlled by the plasma interaction with Jupiter's 
magnetosphere. Europa's O I 1304 A and O I 1356 A emissions are generated by electron 
impact dissociation and are thus a strong function of electron density and temperature. 
Europa's leading side is the downstream side of the plasma convection past Europa. Since 
Europa absorbs the plasma flow including the electrons on its upstream side a rarefied 
wake of electrons is generated on the downstream, i.e. on the leading side. Additionally, the 
electron temperatures are decreased compared to the upstream side due to inelastic collisions 
with Europa's atmosphere. These effect however rather predicts smaller emission near the 
90° west-longitude instead of the observed surplus, but it naturally explains the observed 
smaller emission from the leading side compared to the trailing disk of Europa. 

Europa's UV emissions are intriguingly puzzling. The limitations caused by the HST/ACS 
red leak problem and its associated data analysis problems require further observations with 
HST/STIS at various orbital longitudes and in addition measurements by a possible future 
Europa orbiter to better understand Europa's atmospheric emission and its hemispheric dif- 
ferences. Also further numerical simulations are needed to investigate how plumes combined 
with the time variable plasma interaction possibly control the temporal evolution of Europa's 
UV emissions. 



A. Appendix: Calibration of prism properties 



For the ACS/SBC prism PRL130, both the throughput for wavelengths larger than 2000 
A (the red leak) and the point spread function are not available through the calibration 
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pipeline, so we constrain them here. We estimate both properties with observations of 16 
CygB, a solar analog star. 16 CygB was observed with ACS/SBC PR130L in Cycle 16 (ID 
11325, PI Kuntschner) to improve the existing sensitivity calibrations. In Figure [T41 we show 
observations of 16 CygB dispersed with PR130L. 16 CygB is not a monochromatic point 
source. In case there were no scattering, the dispersed point source would thus lead to a 
'line of emission' on the detector along the dispersion without spread perpendicular to it 
(i.e. in our case without spread in the vertical direction in Figure [H]). The dispersion and 
thus the 'line of emission' is oriented predominantly in the horizontal direction in Figure HH 
(similar to our Europa observations). The scattering additionally spreads the emission also 
to perpendicular to the 'line of emission'. The horizontal direction is also called x-direction 
and vertical direction the y-direction. 



A.l. Long wavelength throughput 



The throughput of the prism PR130L for wavelengths larger than 2000 A can be esti 



mated with the 16 CygB o bservations. We employ a solar analog spectrum (IWoods et al. 



19961 ; iNeckel fe Labsl 1 19941 ; lArvesen et al.l Il969l ; IColina et al.l Il996[ ) which we disperse as a 
point source using the dispersion described with formula ([2]) . To retrieve an estimate for the 
throughput, we integrate the 16 CygB observations along the y-direction on the detector, i.e. 
collapse the 2D emission onto the 'line of emission' and divide it by the modeled emission of 
the solar analog star. The throughput is displayed as a function of wavelength in Figure [15] 
as the red curve. It agrees fairly well with the calibrated throughput of the prism PR130L 
for wavelengths shorter than 2000 A shown in blue. For wavelengths larger than 2000 A, 
our derived throughput serves as basis to construct a model throughput. For A > 2000 A 
we therefore simply use an exponential law which starts at 2000 A with a value of a factor 
of 20 smaller than the calibrated throughput value at A just lower than 2000 A. It reads 



t(A) = to exp 



A - A f 



H, 



red 



(Al) 



with t = 3xl0~ 5 , A = 2000 A and H red = 2000/3 A. Our model throughput for A > 
2000 A (blue curve) matches the throughput estimated from the 16 CygB observations (red 
curve) fairly well. In Figure 



5] we also show for comparison the total ACS/SBC through- 



put rederived by iBoffi et al.l f|2008T) for Sy nphot. In summary, for A < 2000 A we use the 



throughput derived by iBohlin et al.l (120001 ) and for A > 2000 A we use as model throughput 
an exponential law (blue curve in Figure IT5"j) . 



We also apply a modified, i.e. alternative, throughput in our analysis to study system- 
atic uncertainties in the calculation of the surface reflected light. Therefore we assume a 
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throughput dip in the wavelengths range 1800 A to 2300 A by factor of 1/3 applied to the 
model throughput derived in the previous paragraph. The alternative throughput curve is 
shown as dashed blue line in Figure [TSJ 



A. 2. Point Spread Function 

We use the 16 CygB observations (ID 11325, PI Kuntschner) also to construct a point 
spread function (psf) to be used for the Europa observation. Observations of 16 CygB with 
prism PR130L would lead to a 'line of emission' on the detector along the dispersion in case 
there is no scattering. Occurring scattering spreads this line. The scattering is quantified in 
Figure [161 where we show the observations of 16 Cyg B integrated parallel to the dispersion 
in the x-direction (columns 268 to 304) as a function of rows (y-direction) as red curve. 
The observations maximize near row y=432. Near y=432 the flux might be described as a 
Gaussian, but further away it does not show obvious exponential decay and rather resembles 
a power law. 

To model the observed scatter shown in Figure [161 we use the solar analog spectrum, 
disperse it with the prism PR130L properties as detailed in (J2J), and then convolve it with 
a point spread function to be determined. The scattering as shown in Figure HH can be 
modeled by a point spread function of the form 

r / r \ 7i ~ K h 
psf(x,j/) = Co 1+ - (A2) 
L \wJ J 

with r = \J x 2 + y 2 . The mathematical expression in (IA2I) has three free parameters. The 
first parameter w describes the width of the structure. Its value is related to the full width 
at half maximum (fwhm) by fwhm « w X 2 1+1 / K when 7 ^> k. The second parameter 7 
describes the shape of the function around its maximum and the third parameter k models 
a power law decay further from the maximum. The coefficient Co is introduced such that the 
two-dimensional integral over the psf renders unity, i.e. J J psf(x,y) dx dy = 1. This psf 
resembles a modified version of the ^-distribution, a well known phase space density function 
in plasma physics. The ^-distribution approaches a Gaussian for large k and describes a 
power law for small k. The psf in (IA2I) turns into an interesting sub-class of point spread 
functions for 7=2, and w = cts/k. In this case ( 1A2I) leads to 



psf(x,y) = c 



1 + 



-k/2 



(A3) 



This psf resembles for finite k a power law with exponent — k for sufficient distances from 
the source. It turns for large k, i.e. in the limit k — > 00 to a Gaussian with variance a 2 . 
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We convolve the psf (x, y) in flA2j) with the line source calculated with the solar analog 
spectrum discussed above: The coefficients in the psf are constrained by the integrated 16 
CygB measurements shown in Figure [HJ We determine the coefficients of the psf quanti- 
tatively by fitting it (black curve) to the integrated fluxes along the dispersion (red curve). 
We find a good fit for k = 2.65, w = 1.17 and 7 = 20. For these parameters the derived 
psf reproduces the 16 CygB observations well, both near the maximum and along the wings 
(see Figure [T6|) . 

This work is based on observations with the NASA/ESA Hubble Space Telescope ob- 
tained at the Space Telescope Science Institute, which is operated by the Association of 
Universities for Research in Astronomy (AURA), Inc., under NASA contract NAS 5-26555. 
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the Johns Hopkins University during his sabbatical stay in spring/summer 2011. L.R. and 
J.S. acknowledge support by the Verbundforschung Astronomie und Astrophysik durch das 
Bundesministerium fuer Wirtschaft und Technologie. P.D.F., D.F.S., and K.D.R. were sup- 
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Fig. 1. — Europa's O I 1356 emission at the trailing hemisphere. I mage is generated from 
the superposition of all STIS observations on 1999 October 5 (after iMcGrath et al.l 12004 ) . 
Celestial North is up and Jupiter to the left. 




Fig. 2.— Five orbits of HST/ACS/SBC observations on 2008-06-29 (see Tabled]). South is 
to the left and Jupiter is down in each image. 
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Fig. 3. — Solar spectrum after IWoods et al.l (1l996[ ) 



-25 - 




Fig. 4. — Sum of all exposures (top) and modelled reflected solar light from disk of Europa 
(bottom). The same color code is applied as in Figure El North is to the right and Jupiter 
is down in each image. The color scale is identical to the color scale of Figure El 
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Fig. 5. — Observed fluxes (red) along trace (x-axis) summed along y-axis calculated with a 
superposition of all exposures. Modelled fluxes of solar light reflected from Europa's surface 
is shown as black solid line. Reflected light calculated with an alternative model is displayed 
as black dashed line (further explanations in text). Vertical dotted lines show area of O1304 
image and vertical dashed lines displays area of O I 1356 image. 
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Fig. 6. — Model albedo for Europa's leading side as used in our analysis. 
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Fig. 7. — Quantitative comparison of observed fluxes (red) along trace (x-axis) summed 
along y-axis and modelled solar reflected light (black). The difference between the observed 
and reflected light is the emission from Europa's atmosphere shown as black thin line with 
error bars in blue. The emission from Europa's atmosphere is multiplied by a factor of five 
for better visibility (the corresponding values are displayed on the right hand side of each 
plot, respectively). Vertical dotted lines show area of O I 1304 A image and vertical dashed 
lines displays area of O I 1356 A image. 
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Fig. 8. — Average emission of Europa's atmosphere as a function of radial distance from 
Europa's center within columns 272 to 326 have been considered. 
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Fig. 9. — Individual exposures after subtraction of reflected solar light on disk of Europa. 
North is to the right and Jupiter is down in each image. 



31 




-10 -5 5 10 

distance from current sheet in degree 

* Ol 1356 + Ol 1304, extended disk (29 June 2008) 

+ Ol 1356 + Ol 1304, disk only (29 June 2008) 

O Ol 1356 + Ol 1304 (5 October 1999) 

Ol 1356 (5 October 1999) 

A Ol 1304 (5 October 1999) 



Fig. 10. — Total fluxes as a function of nominal distance with respect to Jupiter's magnetic 
dipole equator. Flux on disk and limb for 2008 June 29 observations (r < 1.3i2^ uropo ). The 
dashed lines represent the total fluxes when the alternative throughput model is used (see 
text). 
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Fig. 11. — Ratio of sub- Jovian flux compared to anti- Jovian flux for 2008 June 29 observa- 
tions. For the fluxes within the disk no asymmetry is visible. 
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Fig. 12. — Integrated flux in the x-direction (north south) within the disk of Europa plotted 
as a function of the y-direction (sub- Jovian/anti- Jovian). Red curve are observations and 
thin black line is expected radiation from a radially symmetric atmosphere including the 
point spread function. The black jagged line next to the red one in the lower right panel is 
calculated with the alternative model (see text). 
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Fig. 13. — The time averag e of the squared tidal strain rate for a thin ice shell (see 
Ojakangas fc Stevenson! Il989l . Fig 1), averaged over -90° to 90° latitude and plotted as a 
function of longitude. Values have been normalized so that the maximum value is 1. 



-35 - 




Fig. 14. — Observations of 16CygB. 
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Fig. 15. — Throughput of the prism PR130L . Blue curve is the calibrated throughput for 
wavelengths shorter than 2000 A from iLarsonl ( 20061 ) and modeled with an exponential law 
for wavelengths larger than 2000 A to fit 16 CygB observations. An alternative model for the 
throughput is used in the paper with a reduced throughput by factor 1/3 within 1800 A and 
2300 A (blue dashed line). Red curve is calculated throughput from 16 CygB observations 
with solar anal og spectrum. We additionally show the total ACS/SBC throughput derived 
for Synphot by lBoffi et al.l (120081 ) in their Figure 3 as thin black line. 
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Fig. 16. — Estimation of the point spread function (psf) in y-direction. Observations of 
16CygB in red and model flux to constrain the psf in black. 
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Table 1: Properties of HST program 11186 executed on 2008-06-29 



Orbit 


Start time 


Exposures 


Duration 51 


a 


Total Flux c 


Total Flux c 




in UT 




in s 


in degree 


in 10~ 5 photons s _1 cm" 2 


in Rayleigh - 


1 


04:13:14 


4 


641 


-9.1 


21 ± 2.1 


132 ± 14 


2 


05:48:16 


4 


641 


-1.9 


34 ± 2.2 


216 ± 14 


3 


07:24:12 


4 


641 


6.1 


36 ± 2.2 


226 ± 14 


4 


09:00:07 


4 


641 


9.6 


29 ± 2.2 


187 ± 14 


5 


10:36:02 


4 


652 


8.0 


32 ± 2.2 


202 ± 14 



a per exposure in each orbit 

^Nominal angle of Europa's position with respect to the magnetic equator. 

c O I 1304 A and O I 1356 A combined. The error bars include statistical errors only. 



